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DYNAMICAL BIFURCATION OF MULTI-FREQUENCY OSCILLATIONS 

IN A EAST-SLOW SYSTEM 

A. M. SAMOILENKO*, 1. O. PARASYURt, AND B. V. REPETAt 


Abstract. We study a dynamical counterpart of bifurcation to invariant torus for a system 
of interconnected fast phase variables and slowly varying parameters. We show that in such a 
system, due to the slow evolution of parameters, there arise transient processes from damping 
oscillations to multi-frequency ones, asymptotically close to motions on the invariant torus. 


1. Introduction 

It was shown in the monograph of Krylov M. M. and Bogoliubov M. M. [1] that non¬ 
conservative pertnrbations of a pair of harmonic oscillators under quite general conditions 
lead to the birth of a local attractor homeomorphic to a 2-dimensional torus in a 4-dimensional 
phase space of such a system. The phenomenon of invariant torus bifurcation due to a stability 
loss of a limit cycle when a couple of complex multiplicators cross the unit circle as parameters 
change was examined in ElE] and became widely known with appearance of publications w 
(see also m)- The mathematical framework which enabled us to achieve strict results in the 
analysis of multidimensional invariant tori bifurcations was developed in [M3]. It is also worth 
highlighting references HMH] among other works in this direction. 

The mentioned results are related to the static bifurcations theory, which considers systems 
of the form x = f{x,u), dependent on time-constant parameters u = {ui, ... ,Um), where 
/(•, •): X —)■ is a sufficiently smooth mapping. When one claims, for instance, that 

while the parameters u change along some curve u = u{s), s G (—1,1), there is a stable k- 

dimensional invariant torus being born in such a system due to a stability loss of an equilibrium, 

it means that when s G (—1, 0), the system 

x = f{x,u{s)) (1) 

has an asymptotically stable equilibrium x^ which becomes unstable for s G (0,1), and at the 
same time there exists a continuous (sufficiently smooth) mapping X(-, •): x (—1,1) 

such that X(-, 0) = x* for all s G (—1, 0], and for each s G (0,1) the image of X(-, s): —)■ 

is an asymptotically orbitally stable invariant toroidal manifold of system ([T]). 

Due to certain reasons, it sometimes makes sense to view the family of systems as the system 
in X 

x = f{x,u), u = 0. (2) 

Then one can interpret the bifurcation of invariant torus as follows. The system has an invariant 
set, such that its intersections with planes If^ := {(x,m): u = m(s)} for s G (0, 1) are invariant 
toroidal manifolds % '■= X(T^,s) x {m(s)}. These manifolds shrink to the point {(x*,m( 0))} 
as s —?• -|-0, and each manifold 7^ is a local attractor of system (|2]) restriction to the plane If^ 
(which is clearly invariant). Thus, when one says that the aforementioned bifurcation consists 
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in a birth of an invariant torus in the phase space, it should actually be perceived as a spatial 
phenomenon rather than a dynamical one. 

With the release of publications [20H22], there began the systematic research of truly dy¬ 
namical bifurcations — the effects connected with qualitative changes in a system’s dynamics 
that develop in time and are induced by the actual slow evolution of parameters as they pass 
certain critical values. One of the most resonant achievements in this direction was connected 
with the delayed loss of stability phenomenon in so-called fast-slow systems of the form 

X = f(x, u,e), u = eg{x, n, e), (3) 

where x = (xi,..., Xd) are fast phase variables, u = {ui,, Um) are slowly varying parameters 
and e is a small static parameter. Here we cannot present the complete review of results in 
dynamical bifurcations theory. Let us just note some of the papers [21H27]. In particular, m 
depicts the connection between the delayed loss of stability and “canard”-solutions theory of 
singularly perturbed systems [28], whereas in [26] a dynamical counterpart of the Andronov 
- Hopf bifurcation was studied. Meanwhile, for authors’ best knowledge, the information on 
dynamical counterparts of invariant tori bifurcations seems to be lacking. 

In this paper we consider a (2n -|- m)-dimensional system ([2]) (d = 2n) under assumption 
that its invariant manifold of slow motions (i. m. s. m.) is given by equation a; = 0, i. e. 
/(O, M, e) = 0, and for the linear system 

X = [fx{0, u, 0) + u, 0)] X, 

(which is the system of the hrst approximation for the phase variables x with respect to the 
i. m. s. m.) one can specify the following three zones in the parameters space: the zone of 
asymptotic stability Vg, indehniteness zone "D*, and the zone of complete instability Vu . At 
the same time, the characteristic equation of the operator /j,(0, u, 0) has purely imaginary roots 
for all u from union of the aforementioned domains. In addition, we assume that the system 
ii = £5'(0,M,0) is convergent and its attractor is some point in Under certain additional 
conditions it will be shown that in an (9(-y/£)-neighborhood of the i. m. s. m. of system (|3|) 
one can observe the following dynamical bifurcation. Firstly, while during some time of order 
0{e~^) the parameters u{t) move inside the zone Vg, the phase components of the corresponding 
solution x{t) exhibit exponentially damping oscillations. Next, after u{t) has passed P* and has 
entered Vu, the oscillations’ amplitude starts to grow, and hnally, as f —)■ -|-oo the corresponding 
trajectory of (Ej) is attracted to the invariant torus, asymptotically approaching some trajectory 
on the latter. 

To analyze the system in an 0(y^)-neighborhood of an i. m. s. m., the scaling transformation 
X I—)■ ^/ex is applied, after what the problem becomes a non-local one. The establishment of 
invariant torus’ existence itself does not give rise to many essential complications and is done 
using the same results of u\m in quite the same manner as in However, it is a much 

harder task to determine the non-local attraction basin of an invariant torus of the system 
obtained by the above scaling. For this we have managed to show that the relative measure of 
the attraction domain is estimated from below with a value of order 1 — 

The present article is organized as follows. In Section |2] we formulate a series of conditions 
on the system under consideration, construct its partial normal form in phase variables x and 
make a transition to polar-type variables. In Section [3] the behavior of solutions of the hrst 
approximation system is examined. Section 0] contains our main result, which is based on 
auxiliary propositions of Sections 0] and O concerning the existence of an invariant torus and 
its attraction properties. 
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2. Construction of System’s Normal Form in Phase Variables 

AND Main Assumptions 


From now on we will reqnire system ([3]) to satisfy the following conditions. 

Cl: The right-hand sides of the system are smooth and bounded. Particularly, /(•,•,•) G 
QOO (^2n X Rrn X ^ X X M ^R”^), where C°°{X^y) 

denotes the space of smooth bounded mappings from domain A to set y with bounded 
derivatives of all orders. 

C2: The system has an i. m. s. m. given by equation x = 0, i. e. f{0,u,e) = 0 for all 
(u, e) G R™ X R. 

C3: For every u G R”^ the operator /,^(0 ,m, 0) has purely imaginary eigenvalues ±iuj{u), 
j = 1,... ,n, such that 


inf Ui(u) > 0, inf luAu) — Uk(u)\ > 0, j,k = l,...,n, j ^ k. 


Then for all natural N >2 and s > 2 we may express system Q as 


N 

i ^ Fk{u, e)x’" Fn,s+i{x, u, e)x, 

k=l 


u = e 


■ N 

Gk{u, e)x^ + Gn+i,s{^, u, e) 

_k=0 


( 4 ) 


Here Fk{u, e)x^ and Gk{u, e)x^ are R^"’- and R™-valued homogeneous forms of degree k oix and 
polynomials of degree s and s — 1 in e respectively. The remainder terms of Taylor’s formula 
F]si^s+i{x,u,e)x (here FAr,s+i(x,u,e) is a 2n x 2n-matrix) and G]^^i^s{x,u,e) for ||x|| -|- |e| —?• 0 
satisfy the order relations 


Fisi,s+i{x,u,e) 


O 



Ca+i,s(t, n, £) 


O 




Without loss of generality we can assume that for a hxed natural s and the corresponding 
sufficiently small Cq > 0 the matrix Fi{u, c) is in its real normal form on the set R”* x (—Cqj ^o) 


J(m, e) := diag 


/ eaiiu, e) 


eai{u,e) j 


f 

{u,e) 

-UJn{u,ey 


y ( 

u,e) 


l\ 


Here each function aj{u,e), ujj{u,e) is a polynomial in e of degree not greater than s — 1 and 
s respectively with smooth coefficients depending on u of class C°° (R”*—)-R), and ujj{u,0) = 
ojj{u). To verify this, the following lemma can be used. 


Lemma 1. Suppose that A{-, •) G (R™ x R—, G{-, ■) G C°° (R”* x R—)-R”^), where 
Rdxd d^Yiotes the space of d x d-matrices with real elements. If for all u G R”^ the matrix 
Aq{u) := A{u, 0) has eigenvalues Xj{u), j = 1,..., d, such that 

|Ai(M) - Aj(m)| > 0 Vqj = l,...,d, i^j, (5) 

then for any natural s there exists a mapping T(-, •) G (R™ x R —with the fol¬ 
lowing properties. 1) The mapping e i—)■ T{u,e) is an -valued polynomial of de¬ 
gree s in e with coefficients of class C°° (R™—. 2) There exists £o > 0 such that 
inf(u,£)GR™x(-eo,£o) |detT(M,£)| > 0, and after the transformation x i—)■ T{u,e)x the system 

X = A{u,e)x, u = eG{u,e) (6) 


B{u, e) + e) 


takes the form 


X = 


X, 


u = eG{u, e), 
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where the matrix B{u,e) = is in its real normal form, and 

Bk{-) e C“ , k = 0,...,s, B{-,-) e C°° (M"* X {-eo,eo)^R‘^^'^) . 

Proof. Since condition ([5]) is met, there exists a mapping To(-) G (M™—such that 
the matrix Bq{u) := Tq"^(m)74o(m)To(m) is in its real normal form. Moreover, there exists a 
constant matrix S', with complex elements in general, such that S~^Bq{u)S is diagonal. Let us 
construct a formal change of variables 

X !-)■ '^e^Tk{u)x 
k>0 

with the coefficients Tfc(-) G C°° (R™-—which transforms system ([HD into 

X ='^^e’^Bk{u)x, u = eG{u,e), 

k>0 

where Bk{-) G C°° (R™ —and Bk{u) commutes with Bq{u) for all k>l. 

To do so, let us introduce the dot product {X, Y) := tr(Xy') in and note that as 

{ZX - XZ, Y) = tr [{ZX - XZ)Y] = ii{YZX - ZYX) = - (X, ZY - YZ ), 

the operator X i—)■ ad^X := ZX — XZ is skew symmetric for all Z G thus, = 

kerad^©imad^ (the sum is orthogonal and ad^-invariant). Then, if for any arbitrary Y G 
we denote its orthogonal projection on ker ad^ by Yq, the equation ad^X = Y — Yq with a hxed 
Z has a unique solution X G imad^. In addition to this, if Z has N different eigenvalues, then 
there exists a non-degenerate matrix S (with complex elements) such that S~^ZS is a diagonal 
matrix. Therefore, X G kerad^ if and only if S~^XS is diagonal. 

Next, let and ^>0 formal expansions of A{u,e) and G{u,e) by e: 

respectively. Equating coefficients of powers of £ in the formal equality 

^>0 j>0 i>0 j>0 i>0 j>0 

which is satisfied by Bk{u) and Tk{u), we obtain relations 

Tq{u)Bq{u) = Aq{u)To{u), 

Tk{u)BQ{u) + TQ{u)Bk{u) = Ao{u)Tk{u) + Ak{u)TQ{u) + Rk{u), fc = 1, 2,..., 

in which every operator Rk{u) is dehned by Ti{u), Bj{u), Ai{u) with indices less than k. If we 
assign Tk{u) := TQ{u)Xk{u) for k > 1 and multiply all of the equalities by Tg"^(n) on the left 
side, we will have 

Bo{u) = Tq^{u)Ao{u)To{u), 

-adBo(„)Xfc(M) = T^^{u)Ak{u)To{u) + Rk{u) - Bk{u), k>l. 

Now it is possible to explicitly determine Xk{u) G imadB(,(„) by replacing Bk{u) with the 
orthogonal projection of matrix Pk{u) := (u)Ak{u)TQ{u) + Rkiu) on ker adBQ(.u). (It follows 

from what was mentioned before that Bk{u) equals the diagonal part of matrix S~^Pk{u)S 
multiplied by S on its left side and by S~^ on its right side). 

Clearly, the desired non-formal transformation is T{u,e) = X]fc=o^ 
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Set Sj G to be an eigenvector of matrix Jo{u) := Fi{u, 0) = J{u, 0) which corresponds to 
the eigenvalue iuj{u), k = 1,... ,n. Since 


Jo{u) = diag 


( ° 


0 


/ 0 -UJn{u)\ 

\ujn{u) 0 yj ’ 


the vectors Sj are independent of u. Let us compose a matrix S', whose first n columns are 
the vectors Si, ..., and last n columns are their complex conjugates respectively, and let us 
define basis forms 

‘^q(2/) := ei,q(l/) = <^q(l/)Si, (7) 

where q := (gi, • • •, q2n) e = xf ■ ■ ■ xf”. 

We can now proceed to the construction of a transformation that converts the W-jet of 
system (jlj) to its normal form in fast variables under an extra assumption of the absence of 
resonances up to a certain order between frequencies u}k{u), k = 1,... ,n. To make the corre¬ 
sponding statement, we dehne an n x 2n-matrix I = [En, —En], where is the n-dimensional 
identity matrix, assign 


u{u) := .. .,Un{u)) , |q| := |gi| H-h \q 2 n\ , 

denote the z-th unit vector of the coordinate space (i. e. the vector, whose z-th coordinate 
is equal to 1 and the rest are zeroes) by e*, and introduce the following sets for positive numbers 
u and a 

Ai{N, I/) := G : |(a;(M), /(q — 6 ^)) | > z/ Vq G : 2 < |q| < W, J(q — e*) 7 ^ O} , 
Ao{N,iy) := {u G M"*: |(a;(zz),/q)| > z/Vq G Z^”: 2 < |q| < iV, Jq ^ O} , 

n 

A(N,,.):=f\At{N,F), <i}, Bf := Bf (0), 

i=0 

9fo(W) := {q G Z^": 0 < |q| < iV, Jq = 0} , iHi(iV) := {q G Z^": 2 < |q| < iV, /(q - e^) = O} . 


Proposition 1 . Suppose that conditions Cl C3 are met, and s >2 is a fixed natural number. 
Also, for some W G N, iV > 2 and u > Q let the set A{N, v) he non-empty. Then there exist 
numbers 5 > 0 and £o > 0 such that after the change of variables 


N 


N 


X = 


z/ + ^Wfc(n,£)/, u = v + e'^Uk{v,e)y'', {y,v,e) e B]''x A{N,u) x {-eo,eo), {8) 


k=2 


k=l 


where 


and 


S —1 


Xk{v,e)y^ = Xk,j{v)y^, U{v,e) = ^e^Ukj{v)y'', 

j=0 j=0 


Xk,j(-)y^eC^(A(N,iA)^m^^), UkA-)y^eC^(A(N,iA)^m^) VyeR^fi k = l,...,s, 
system (jl]) takes the form 


2n 

y = J(n,e)z/ + E E i7i,q(n, £)ei,q(z/) -h HN,s+i{y, v, e)y, 

i=l qG9Ifc(V) 


V = e 


^ <^q(z/)Cq(z;, e) + C 7 v+i,.(z/, v, e) 

qe9Io(N) 


( 9 ) 
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Additionally, 


s-l 


C^{v,e) = 

j=0 j=0 

where ^ C°° (^(iV, z/) —)■ C), Cqj(-) G C°° {A{N , u) ^C^), and the remainder terms 

of Taylor’s formula H]^^s^i{y,v,e)y and C7v+i,s(i/, u, e) satisfy the order relations 


HM,s+i{y,v,e) = 0 [\\y\\ +e 


\N . s+1 


CN+i,s{y,v,e) — O (+ £"") , ||i/|| + |e:| ^ 0. 


Proof. Let us apply the polynomial transformation ([8]) to the truncated system 


N 


N 


X = J{u, e)x + Fk{u, e)x^, ii = e Gk{u, e)x^, 


k=2 


k=0 


If the resulting system is 

y = J{v, e)y + ^ Hk{v, e)/, v = e^Ck{v, e)y^, 


k>2 


k>0 


then the following equalities must hold 


N 


j(v, £)y + Hk{v, e)y^ + ^ 

i=2 


d{Xj{v,e)y^) 


k>2 


dy 


J{v,e)y + '^Hi{v,e)y^ 


i>2 


+ 


N 


d{Xj{v,e)y^) 


dv 


+ 5 : 

i=2 

= J iv + e^Uj{v,e)y\ 

V j=i 

N / N 

lv + e'^Ui{v,e)y^, 


i>0 


N 


+ 


i=2 


Z=1 


N 


J 2 Ck{v,e)y^ + J 2 


k>0 


i=i 


d[Uj{v,e)y^] 
dy 


'^Ci{v,e)y 
!0 

N 

y + ^Xi{y,e)y 

i=2 

N 

y + '^^i{v,e)y 
i=2 

J{v,e)y+ ^Hi{v,e)y 


i>2 


+ 


N 

E 

i=i 


d{Uj{v,e)y^) 

dv 


^Ci{v,e)y^ 


i>0 


N 


N 


'Y^Gj (u+gy^ Ui{y,e)y\i 

j=0 


1=1 


N 


y + "^^i{'vx)y" 


i=2 


Equating the homogeneous forms of y in the left- and right-hand sides, we obtain 

C'o(u,e:) := Go{v,e), 

Gi{v, e)y -f Ui{v, e)J{v, e)y + ^^^ Gq{v, e) = 

^ ^ . d (Xk(v,s)y^) ^ ^ d (Xk(v,£)y^) 

ffk(v,£)y -I - ^ - hJ(v,£)y + £ - — - -Go(v,£) = 

= J(v, £)Xk(v, £)y^ + Fk(v, £)y^ + Mk(v, £)y^, k = 2,...,N, 
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CkKy,e)y H--- J[v,e)y^e ---C'o(t,£) = 


dy 


dv 


= £)/ + Gi{v, e)Xk{v, £)/ + Gkiv, e)y^ + Nk{v, e)y^, fc = 2,..., iV, 

ov 

where all of the forms Mk{u,e)y^ and Nk{u,e) are determined by the forms contained in the 
original system, the resulting one and in the transformations and have indices less than k. The 
forms within the given equalities can be expanded as Hk{u,e) ("NJ j{u) and so on. 

Having introduced operators 

d ■ d ■ 

^Jo{v)y ■ ■) ■ .= -^^J^itj^y 


and having equated coefficients in left- and right-hand sides, we arrive at the homological 
equations for determination of Xkj{u), Hkj{u), Uky{u), Gky{u): 


Uij{v)Jo{v)y = Gij{v)y - ^ Ui,i{v)Jj-i{v)y+ 

i=0 

where j = 0,..., s — 1, and 

^Jo{v)yXk,j{v)y^ = Fkj{v)y'' + Pk,j{v)y'' - Hkj{v)y'', 

djo{v)yUk,j{v)y'' = Gi,o(u)Xfcj(n)/ -F Gk,j{v)y^ + Qk,j{v)y'^ - Gk,j{v)y^, 

where k = 2, ... ,N, j = 0,..., s. Here the forms Pkj{u), Qk,j{u) are constructed using the 
forms found from the analogous equations on the previous step, i. e. from the equations in 
which there is an index k — 1 instead of k. The forms Hkj{v) and Gkj{v) are chosen in such 
a way that the resulting system has, in some sense, as a simple structure as possible. Since 
Jo(n) is non-degenerate, Uij{v) is explicitly found by assigning Gij{v) = 0, j > 0. After that, 
we switch over to determining other required forms by replacing them with their expansions in 
the basis forms (]7|): 


dGoj-i-i{v) dUi,i{v)y 

Ui^i{v)y -- Goj-i-i{v) 


dv 


dv 


-Gij{v)y, 


2n 


^k,j{v) — Ari,qj(n)ei,q(i/), Uk,j{v) = ^ ^q(i/)f/qj(n), etc. 

|q|=fc *=1 |q|=fc |q|=fc 

As S~^Jo{v)S := diag [ia;i(n),..., ia;„(n), —ici;i(n),..., —ia;„(n)], it is not hard to deduce the 
equalities 

9Mv)y^^iy) = [^^iy)t Mv)y = ^ 


dt 


t=o 


(S = i(u;(v),Iq)^q(y), 

^Jo(v)yei,q(y) = [ei,q(2/)]' Jo(v)y - Jo(v)ei^q(y) = ^ = 


= i(a;(n),/(q-ei)ei,q(i/). 

But then the equations for determining the desired forms’ coefficients become 
i -f(q ^0) 

i (a;(n), Jq) Uqj{v) = Gi,o(n)Xqj(n) + Gqj(n) + Qqj(n) - C'qj(n). 

If n G Ai{N, z/), i 7 ^ 0, then in case /(q —Sj) = 0, we can declare Hi^q^iy) = Fj qj(n) — Pj^qj(n), 
~ otherwise Hi^qj{v) = 0, and at the same time we will hnd Xjqj(n). Similarly, 
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if u e ^o(^, I^), then C'qj(n) = Gi,o(n)Xqj(n) + G'qj(n) + Qqj(n); f/qj(n) = 0 if Jq = 0, and 
Cqj(n) = 0 if Jq 7 ^ 0, in which case we explicitly hnd f/qj(n). 

Having performed the constructed change of variables in (jl]), we obtain system ([Qj). □ 


Let us now switch over to complex variables z = {zi,, Zn) € C"" in system ([9j) with the 
substitution 


n n 

y = Y^ = 2Re 

i=i i=i 


n 

-i=i 


System ([9]) takes the form 


= 


eaj{v,e)+ iu:j{v,e) + ^ hj^^{v,e){\z\f'^ 

3<2|p|+l<Af 

+0 \\z\\) , j = 1,.. .,n, 


Z4-\- 


V = e 


^ c^{v,e){\z\ 

0<2|p|<V 


|2P + 0 ^ ll~ll^+^ 




where p e Z" and hj^^{v,e) := (p,p)+e^.(n, e), Cp(n,£) := C(p,p)(n,£), (p, p) : = 

(pi,... ... ,p„), (|z|) = (|zi|,..., |z„|). The order relations O + llpllj 

and O denote the remainder terms of the same kind as Hf^^s+i{y-,v.,e)y and 

<5a+i,s(i/, u, e) respectively. One can also easily ensure that the equations for Zj are complex 
conjugate with the equations for Zj. 

Throughout the rest of this paper we will assume that the parameter e is non-negative. 
Having introduced the polar-type coordinates r^, ipj\ mod 27r by Zj = , j = 1,... ,n, 

having dehned aj-p(n,e:) := Rehj-p(n,£), 6j-p(n,e:) := Imhj-p(n,£), r = (ri,...,rn), y/r = 
{y/rl ,..., (p = (pi,..., pn) and having assigned s = {N + l)/2, we come to the system 


O = 2^ 


aj[v 


,e)+ ^ e'Pl ^aj'p(n,e)rr 

3<2|p|+l<N 




V = e 


£|plcp(n, £)rP -1- n, <p, e) 

0<2|p|<A 

= Qj{v, £) + ^ £'^'^i,p(n, e)r^ + n, p, e), j = 1 ,..., n, 

3<2|p|+l<V 

where the remainder terms can be written as 

Rj{r,v,(p,e) := ^ Oj-q(v^, n, p, ejxA’'* + 2 ^ aj,q(v^, n, p, e) 

|q|=A+l |q| = l 

^j{r,v,ip,e) := ^ 6 ^-q(^/if, n, p, ^ 6 ^-q(^/if, n, p, 

|q|=JV+l |q| = l 

Z(r, n, p, e) := ^ Cq(v^, n, p, + c(x/if, n, p, e). 

|q|=V+l 


( 10 ) 
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Here the functions aj,q(p, n, e), aj,q(p, n, e),... are smooth in [0, ^o]" x K™' x T" x [Oj^o], 
and > 0, £o "C 1 are some positive constants. 

If we declare for a couple of vectors p = (pi,... ,Pn), <i = (?i, ■ ■ ■, ?n) an operation p • q = 
{piqi,.. .,Pnqn) and assign 

a{v) := (ai(n),... ,a„(n)), A{v) := - , c{v) := Co(n,0) 

(here Sj denotes the j-th coordinate unit vector of the space M"'), it will allow us to rewrite 
system 00 as 

r = 2e [a(n) — A{v)r + eB{r, v, e)] • r + • R{r, v, <p, e), 

V = ec{y) + e^W(r, v, e) + v, (p, e), (11) 

ip = u:{v) + e'^{r,v,e) + • ^{r,v,p,e). 

All of the functions which appear in this system are bounded in [0, p]” x MB x T” x [0, ^o], where 
0 < g < qo/eo, and their smoothness properties are determined by the corresponding terms 
of system flT0|) . Furthermore, the order relations \\R{r,v,p,e)\\ = Odlx/Lll), ||<F(r, n, <p, e)|| = 
0(||i/r||) hold uniformly in n G p G T", e G [0,£o] when ||r|| —)■ 0. 

For the sake of simplicity, we shall consider the case when the domains Bg, P* and which 
were mentioned in the Introduction, are formed by nested balls. More precisely, let us introduce 
the following notations for a triplet of numbers Rq,R^,, R* such that 0 < Rq < R^ < R*: 

ao := min inf aj{v), a* := — max sup ctjiv), a* := sup ||a(n)||, 

min (kl(n)e,0 , sup ||A(n)|| 

and state some additional assumptions. 

C4: There exist such numbers Rq < R^ < R* that a* > 0, cio > 0, A* > 0. 

C5: The conditions of resonances absence are fulhlled: B'^t C A.(iV, u) for some N > 3, 
u > 0. Besides, if < 5, then 0 G A,(5, z/). 

C6: The conditions of system v = c(v) convergence are met: there exists x > 0 such that 
(c(v),v} < —x||n||^ for all v G B^,. 

C7: All components of vector r* := A“^(0 )q;( 0) are positive. 

Furthermore, without loss of generality we can assert that r* = (1,1,..., 1). This can always 
be achieved using the scaling transformation r >—> r* • r. 

In accordance with condition C4, we now take that 

Vg = B^.\BZ, V. = BZ\B-, Vn = B-. (12) 

3. Analysis of the First Approximation System 

In order to have at least rough understanding of how system dm solutions behave, let us 
focus on the hrst approximation system 

r = 2e[a{y) — A{v)r] • r, v = ec{v), p = oj{y) + e^{r,v,e). 

Currently we are mainly interested in dynamics of the subsystem for the variables r, v 

r = 2e[a{v) — A{v)r] • r, 

V = ec{v). 

Proposition 2. The origin is a global attractor of system ms within a ball Bf}*. 


(13) 

(14) 
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Proof. Condition C6 implies that along every solution of system flT^ the function {v, v) tends 
to zero monotonously. □ 

Let {T(f)}^>g be a forward trajectory of system flT^ and let us consider the solution r(-) of 
the system 

r = 2e[a{v(t)) — A{v(t))r] • r, 

such that r(0) G (0, cxd)"'. Since the j-th component rj{-) of this solution can be viewed as a 
non-trivial solution of a linear homogeneous equation with a continuous coefficient, it yields 
rj{t) > 0 for all t >0. 

Proposition 3. Suppose that n(0) G Ttg o,nd T* = supjf > 0: v(t) G Ps}- Then |r(f)| < 
|r(0)|e-2="** for all t e [0,T*]. 

Proof. By condition C4 on the interval [0, T*] we get 

< 2e [-a* |r(f)| - (y4(n(f))r(f),r(f))] < 2e [-a* |r(f)| - ||r(f)||^] < -2ea^ \r{t )\, 

which provides the estimate for |r(f)|. □ 

Corollary 1. For all v &T>s the derivative of |r| along trajectories of subsystem fll3p does not 
exceed —2ea^ |r|. 

Proposition 4. The point (r*,0) is a global attractor of system in the domain 

(0,cx))" X 5^*. 

Proof. It follows from the inequality 

< 2e [\a{v{t)) • r{t)\- {A{v{t))r{t),r{t))] < 2£||r(f)|| [a* - A, ||r(f)||] 

that |r(f)| is decreasing while ||r(f)|| > a*/A^,, and hence, it does so for at least as long as all 
points of the hyperplane |r| = |r(f)| stay outside the sphere ||r|| = a*/A^. Or in other words, 
it is decreasing while the distance from the hyperplane |r| = |r(f)| to the origin is greater than 
a*/A^. Since this distance equals |r(f)| /^/n, then no matter how small 5 > 0 we choose, there 
will be a unique positive moment of time starting from which r{t) belongs to the bounded set 
{rG(0,cx))”: \r\ < ^/n {a* + 6) /A^}. 

Meanwhile, as soon as at some moment 0^0 the point v{t) enters Pu, the inequality 

> 2e [a* |r(f)| - A* ||r(f)||^] > 2e \r{t) \ [a* - A* \r{t)\]. 

will become valid. Consequently, starting from some moment of time ti > to the inequality 
|r(f)| > (a* — 5) /A* holds. Thus, if we declare 

lC-.= {reWf-. (a, - 5)/A* < |r| < [a* + 5) /A,] , (15) 

there will be a moment of time t/c = tic{r{0)) > ti such that G /C, and then r{f) ^ K, for 

every t > tic- It is noticeable that by choosing the number Sq small enough, without loss of 
generality, we can thereby suppose that g > y/n{a* + 5) /A*, and accordingly, |r| < g for all 
r G /C. 

Let us prove that r{t) —>■ r* when t —)■ +cxo. Consider the limit system 

r = 2e[a{t)) — 74(0)r] • r. 

It has a positive definite Lyapunov function in (0, cxo)” (relative to the equilibrium r*) 

n n 

:= X] - 1 - Inc) = |l| - ln]^ri - n, 

2 = 1 2 


(16) 
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with a negative definite derivative along the limit system. Indeed, 

(VVo(r), 2£[a(0) - A(0)r] • r) = 2^ ^ ~ "4(0)r]. = 

= 2e {r — r*, a(0) — 74(0)r) = —2e (r — r*, ^4(0) [r — r*]) < —2ey4* \\r — r*\f . 

Now let us assign q := supo<||^||</j* ||u||~^ [||ci('^^) “ Q!(0)ll + ll^('^) “ ^(0)11 ^'] compute the 
derivative along trajectories of system of the function 

l/(r,n) :=l/o(r) + A||n||V2, (17) 

where A > / {2A^,>c). Sylvester’s criterion of positive definiteness of a quadratic form claims 

the existence of such a number /i > 0 that 

(Vho(r), 2e[a{v) — A{v)r] t) + e\ {c{v),v) = —2e {r — r*, 21(0) [r — r*]) + 

+2e {r — r*, a{v{t)) — a(0) + [2l(0) — y4(n(f))] r) + eA {c{v),v) < 

< —e [2A^ ||r — r*\f — 2q\\r — r*|| ||n|| + Ax ||n||^] < —e/i [||r — r*\f + (18) 

for all V G and r G (0, cxd)"' such that |r| < g. Since {r{t),v{t)) G /C x for all sufficiently 
large t, then V{r{t),v{t)) — )■ 0 for f —)■ +cxd. But thus r{t) —)■ r* when t —)■ +cxd. □ 

Corollary 2. For all v G and r G (0, cxd)” such that |r| > \/n{a* + 5)/A^ the derivative of 
funetion |r| along trajectories of subsystem ffldj) does not exceed —2e6 {a* + 5) jA^. If, otherwise, 
V G T>u and 0 < |r| < (a* — 5)/A*, this derivative is greater than 2e6 |r|. The set K, x is 
a forward invariant set of system dH-dH. Moreover, each forward trajectory of this system 
such that (r(0),n(0)) G (0, oo)"' x enters K. x 

Let J(r) denote the Jacobi matrix ^ [(o(0) — A(O)r) • r] and Hy^{r) be the Hesse matrix of 
the function Vo(-) at the point r. One can easily verify that the quadratic form {Hy^{r*)r,r) is 
positive definite. 


Proposition 5. The linear system r = J{r*)r is asymptotically stable and the derivative of the 
quadratie form {Hy^{r*)r,r) along trajectories of this system is negative definite. 


Proof. We have the inequality 

— (r — r*, H(0) [r — r*]) = (VV))(l), [a(0) — y4(0)r] • r) = 
= {Hyoir*)ir - r*),J{r*){r - r*) + 0(||r - r*||^). 


Since the left-hand side is a quadratic form, the right-hand side should also be one. 
concludes that 


{Hy^{r*)r, J{r*)r) = — (H(O)r, r) Vr G M"". 


This 


□ 


4. The Main Theorem 

Everywhere in what follows, we assume that g > max{-yn {a* -|- J) /H*, 1}, and therefore, 
the set /C given by formula flT^ is contained in the simplex 

Sgy := {r G M+: |r| < p} C [0, gf . 

Before we proceed to the formulation of the main theorem, let us point out an important 
property of a union of sub-level sets of the function Vo(‘); defined by 0161) : 

|Jv'„-‘(|o,c]) = (o,«d)”. 

c>0 
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This fact, particularly, is a consequence of the following simple lemma. 
Lemma 2. For k > 0 and c > 0 consider the set 


Qeik, c) := {r e M”: > e j = 1,..., n} . 

If Eq G (0,1) and c> 1, then 

([O, lineal +c- l]) := {r G (0,oo)"': Vo{r) < |hie^| +c- l} C Qs{k,c) We G (0,£o)- 


If additionally 


hie 


k/n 


> p — In p, then 


Qe 



ns,crp([o, [in£‘|]) 


Ve G (0,Co)- 


Corollary 3. The estimate mes (5^ \ Vq ^ ([O, |lne^|])) = holds when e —)■ +0. 

Now we can proceed to the statement of our main result. 


Theorem 1. Suppose that conditions C4 C7 are fulfilled and 0 < k < N — 2. Then there 
exists such So > 0 that: 1) for every e G (0,eo) th^ solution {r{t),v(t),ip(t)) of system flTT]) with 
the initial condition (r(0), n(0), <^(0)) G x x T” can be extended to the semi-axis [0, cxo), it 
satisfies the inequality |r(f)| < |r(0)|e“^"** on the interval [0,Ti(e)) := {t > 0: v{t) G T)s}, and 
there exists such an instant T 2 (e) > Ti(e) that r(t) G 1C, v(t) G T>u fort > T 2 {e), where sets 1C, 
T>s and T>u are defined by flTSl) . (IT^ .- 2) system flTT]) has an n-dimensional invariant torus % 
located in an 0{e)-neighborhood of the torus {r*} x {0} x T", and the system’s restriction to 
Ts has the form (p = Ci;(0) + ef(95, e) where f(-, £):¥"—>■ M" is a Lipschitz vector field; 3 ) if in 
addition r(0) satisfies Vo(r(0)) < |lne^|, then there is a trajectory {{fit),v{t),p{t))}^^,^ on the 
torus Te such that 


lim [|r(f) - f(f)| + |n(f) - h(f)| + |v?(f) - p{t)\] = 0. 

t—>-+oo 

The statement remains valid for arbitrary r(0) G Sg if the additional conditions hold: 
Rj{r,v,:p,£)l_^= 0 y{r,v,ip,£) e Sg x 5^. x x (0,eo), j = 

The proof of this theorem is drawn from the statements of Sections |5] and O 


(19) 


5. Preliminary Analysis of the Normalized System 

In what follows assume that conditions of the main theorem are met. The following propo¬ 
sition captures a series of similarities between the dynamic of the first approximation system 
and the behavior of system flTT]) . 

Proposition 6. Let {r{t),v{t),(p{t)), t ^ I, be an non-extendible solution of system flTT]) such 
that r(0) G Sg, n(0) G Bff. Then for sufficiently small £o > 0 and every £ G (0,£o) Ihis 
solution has the following properties. 1) The interval I contains the positive semi-axis, and 
hence, SgXBfj, xT” is a forward invariant set of system flTT]) . 2) There exists a moment of time 
A > 0 after which v{t) does not leave some 0{£)-neighborhood of the origin o/R™'. Moreover, 
on the interval [0,re] the function ||n(f)|| is monotonously decreasing. 3 ) While v(f) G T>s, the 
function |r(f)| decreases with an exponential rate and it satisfies inequality |r(f)| < |r(0)| 

4) The set K, x T>u x is a forward invariant set of system (Hj, and there is a non-negative 
moment of time, since which {r{t),v{t),(p{t)) & 1C x T>u x ¥"■. 
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Proof. By computing and estimating using Corollaries [H [2] the derivatives of the functions 
|r| and {v,v) along trajectories of the corresponding subsystems of system ca. it is easy to 
verify that under the condition that Eq is sufficiently small in the corresponding domains these 
derivatives have the same signs as the derivatives of the functions |r| and (u, v) along trajectories 
of system flT^ - (IT^ . In the same fashion as in proofs of Propositions [21 [3l HI we obtain the 
desired result. □ 


The presence of the term • R{r,v,^p,£) in system (ITTl) makes it harder to establish 

a counterpart for Proposition jH The next proposition provides restrictions on initial values of 
system flTT|) solution which guarantee that, starting from some moment of time, this solution 
enters and remains inside an 0(\/£)-neighborhood of the torus defined in the phase space by 
equations r = r*, v = 0. 

Proposition 7. There exist such positive numbers Eq and C* that for all e G (0, £o) 
solution (r(t),v{t), ip{t)) of system (ITT]) with the initial values r(0) G 5^ fl 1/^“^ ([0, |ln£^|]) and 
u(0) G Bfft there is a moment 4 > such that 

||r(f) — r*\\ < C*y^, ||u(f)|| < C*E Vf > te- (20) 

If additionally conditions (unD are met, then the existence of te is guaranteed for any solution 
of system (ITT]) such that |r(0)| < g, rj{0) >0 (j = 1,... ,n), u(0) G 

Proof. Turning back to Propositional it is enough to verify the first inequality (1201) . Let Cq 
be a constant that bounds from above each of norms ||i?(r, u, e)||, \\R{r,v,ip,E)\\, ||TP(r,u,£)||, 
||Z(r, V, if, £)|| on the set [0, pj” x Bfl* x x [0, £o]- By (IT8]) on the set (0, p]"' x T"- x Bf^, x [0, £o], 
the derivative of the function V (r, v) (see (IT7]) j along trajectories of system (ITT]) can be estimated 
as 

■= (VV'o(?’), 2e: [a(u) - A{v)r + EB{r,v,E)] •r + • R{r,v, p,e)) + 

+A (u, Ec{v) + E^W{r, V, e) + v, (p, e)) < 

< —EH (||r — + ||u||^) + 2e^ (r — r*, B{r, v,e)) + e^^"^ (r — r*, • R{r, v, (p, e:)) + 

{v, W{r, V, e)) + {v, Z{r, v, <p, e)) < 

< —E [fi (||r — r*\f + ||u||^) — ||r — r*\\ [2 CoE + ^)||) ~ 2£ACo ||u||] . 

Let us show that for sufficiently small Eq this derivative does not exceed a certain negative value 
on the set 

|(r,u, 9 ?,£) G [Qe{k,c) n5^] X Bff, x T” x (0,£o]: \J\\r - r*\\^ + ||n||^ > OV^Co/pj , (21) 


where c > A[i?*]^/2 + 1. First of all, one can notice that 

||u|| [h ||u|| - 2eAC'o] < e^X^CI/h Vu G Bff,. 

Let {r,v,ip,E) belong to set (12T]) . If ||r — r*|| > 1/4, then inequalities Vj > e~^E^ yield that for 
the hrst n coordinates of the point of set (l2T]) we get 

^A/2-i ||^-i/2 . u, (p, £) II < e^P-^Coe^pE-’^P = £(^-^-2)/2Coe^/2, 


and then for sufficiently small Eq 


V^{r,v,p,E) < -- 


^ - 2CoE - £(^-^-2)/2C'oe"/2 


+ £ A Cg //i < 0 G (0, £o] • 
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If, on the contrary, ||r —r*|| < 1/4, then |rj — 1| < 1/4, and thus, Vj > 1/4. Taking into 
account that in this case • i?(r, n, (^, e:) || < 2Cq, for sufficiently small £o and all e G (0, eo] 

we obtain 

^ [^^ (II’’ “ ’’11^ + ll^lfl “ 4 V£C'o ||r - r*|| - 2eXCo ||n||] < 


< -£ 


/i 



< 0 . 


In accordance with Lemma |5] the set Vq ^ ( [O, [ine’'! + c — l]) lies in Qe(A;,c), which means 
that 

© ;= y-' ([0, line'll + c - 1]) n [5, x ] C [Q, {k, c) n x , 

and Proposition [6] together with the estimates for V^{r, v, cp, e) imply that the set © x T” is 
forward invariant. Furthermore, if (r(0),n(0)) G ©, then there exists a moment r* > 0 such 
that V{r{t)^v{t)) < c*{e) for all t > t*, where 



In fact, in a closed ball centered at (r*, 0), the function !/(•, •) reaches its maximal values only 
on the boundary. Therefore, in an open ball of radius 6\/^Co/h) which is centered at (r*,0), 
the function !/(•, •) takes values that are less than c*{e), which means that this ball lies inside 
V~^ ([0,c*(£)]). Now the existence of 4 results from the negativity of V{r,v,ip,e) at points of 
set ([21]) • 

Obviously, if r G 5^, Vo(r) < |ln£^| and ||n|| < R*, then V{r,v) < [ine’'! +c—1. Hence, the 
set © contains the set ([O, |lne^|]) O 5J x 

Furthermore, since Vo(r) ~ | {Hy{r*){r — r*),r — r*) near r*, one can specify such C* > 0, 
that the set H“’([0, c*(e:)]) lies in a ball of radius C*^ with center at (r*, 0) for all sufficiently 
small £ > 0. It enables us to assign = max{r£, r/}. 

Finally, if conditions flT^ are met, then Rj{r,v,(p,e) = y/r^Rj{r,v,(p,e), j = and 

the proof runs as before in case r(0) G 0 (0, cxo)” if subsystem for r (ITT]) is replaced by 


r = 2e 


a{v) — A{v)r + eB{r, v, e) + ^R{r, n, (p, e) 


• r. 


□ 


Remark 1. It follows from Propositions [6] and [7| that when the forward trajectory 
[J^>o(r(f), n(f), <p(f)) has no common points with the set © x T”, then there is a moment 
of time after which r{t) G /C \ Vq~^ ([0, |ln£^|]). 

Hereafter, our main question will be whether system Dll possesses an invariant torus close 
to the invariant torus of the first approximation system and if so, what its basin of attraction 

is. 


6. Existence of an Invariant Torus and its Basin of Attraction 

Taking into account the already proved propositions, we will conduct further analysis of 
system ffTTj) in a neighborhood of the torus r = r*, v = 0. 

In order to simplify our notations, we will combine the variables r and v into one vector 
variable y = (r, v) (this (n+m)-dimensional variable has no connection with the 2n-dimensional 
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local variable in Section |2]). Let us rewrite system cm as 

y =eF{y,e) + e^/’^G{y,Lp,e), 

0 =^{yG) + £^^'^H{y,^,e), 

where 

F{y, e) := (2 [a(n) — A{v)r + eB{r, n, e)] • r, c{v) + eW (r, v, e )), 

G{y, := (v^ • R{r, V, (f, e),e^l^Z{T, n, <^9, e)) , 
bj{y,e) := a;(n) + £^(r, n, e), H{y,ip,e) := • <h(r, n, v?, e), 

and let us assign y* := (r*, 0). Then F{y*, 0) = 0. Since 

0 c'(0) ) 

and according to Proposition [3] and condition C6 the linear systems r = J{r*)r and v = c'(0)n 
are asymptotically stable, which means that all of the eigenvalues of these systems’ matrices 
have negative real parts. Consequently, the system y = F'{y*, 0)y shares this property, too. It 
is well known that there exists a positive definite quadratic form which has a negative definite 
derivative along trajectories of an asymptotically stable linear system with a constant matrix. 
This positive definite quadratic form sets a dot product structure. Therefore, we will further 
assert that the space is endowed with the dot product (•, •) for which the quadratic form 

(^Fy{y*,0)y,y'^ is negative definite. 

Now, we can choose positive numbers 7 , a and £0 in such a way that the inequality 

( [Fy{y, e) + (p, e)] z, z) < -27 (23) 

V(i/,^,V?,£) e X X T” X [0,£o]. 

holds. It implies that x T” is a forward invariant set of system (l22l) . Thus, for each 

point (i/, (p) G B^~^"^{y*) x T” its forward trajectory, denoted by {{ritiy,^)),4>t{y,^)}t>oy lies in 
B^~^^{y*) X T”. In other words, in B'^~^^(y*) x ¥"■ system (1221) generates the semi-flow 

-ST"!!/") XT" XT” 

It should be mentioned, that actually, as it follows from results of Section [31 system (122)) 
generates a semi-flow on the set Sg x B^l* x Moreover, each point that enters the set 
[5«nV'7(|0,lln£*|l)] X B^. X T” at some moment of time under action of this semi-flow 
necessarily enters B^~^^{y*) x T” after some instant of time and later on keeps moving inside 
0 (\/£)-neighborhood of the torus {y*} x 

Now, let us show that for any sufficiently small e > 0 the set B^^'^{y*) x contains an 
n-dimensional invariant torus of system ( 122 )) which attracts all forward trajectories of this set, 
and hence, the attraction basin of this torus contains the set [iSg n Vq~^ ([ 0 , k |ln£|])] x B^^ x T”. 

Remark 2. If a is small enough, condition C5 yields that B^ C M(5, v) and then for v G S™ 
results of Section [2] remain valid for N = 5. Because of this, we will further consider system ( 1221 ) 
with N > 5. 

It follow from (123)) that the evolution matrix hi* of the linear system 

i = e [Fy{yt{y, ‘f),e) + e^^~^'>/^Gy{rit{y, (p), My, P),^)] ^ =■ y, P, 
can be estimated as 

< g-2£7-(i-*), t>5>o. 
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One can choose the positive nnmbers a and K in such a way, that for all (?/*,</?*,£) G 
X T” X [0,£o] fhe following inequalities will hold 


p -0j{y\E) - (y , e) II < 


V 




(24) 


<A'[||»'-»=||+e«/= 

||F(!,‘, e) + ,f\e)- F{y\e) - £)- 

- K(y,E)(»■ - »')ll < 

- [Py(y^, c) + £'^“^*''^G;,(!/^ ip^, £)] (s' - y^) [I + 

+ II £) - U(!/t s) + £'"-"’'"G;(s^ p\e)- £'''-'IGg;(s', £)] (s‘ - s")|| + 

^j(W-2)/2 ||g(,,i_^i_ 5) -G(s',s>t£)|| < 

< K [(lls' - s^ll + ||s'' - s“|| + IP^ - V2“||) p' - y|| + \y - S>"||] . (25) 

Proposition 8. Let us declare := {{y-,z) G x y^z ^ ||i/ — z\\ < e} 

and M :=4i^/7. There exists such Eq > 0, that when (0,£o) to each point {y, z,(p) gS^xT” 
there is a unique corresponding point 6{y, z, cp) G such that for all t > 0 the inequalities 

llht {y, t) - Vt (z, d(y, z, <^))|| < 2e-^^* \\y - z||, 


Ut {y, t) - <Pt {z, e{y, z, <^))|| < —e \\y - z 


(26) 


hold. Furthermore, G C{Bs x T”^T”); for every {y,z\p) E B^ x T”, i = 1,2, it is 

true that 

11^ (l/, t)-G ( 2/2 V?) II < y 11^^ - > 

and for any fixed point {y, z) G Be the mapping 6{y, z, •):¥"—)• T" is a homeomorphism. 

Proof. Set Tie fo be the space of continuous mappings 

Ue := R+xBeXT-3 (f, y, z, p) ^ (C(t, y, z, p),'iP{t, y, z, p)) G x T" 

which for all (t, y, z, p), {t, y, z\ p) G M+ x Be x T”, z = 1, 2 satisfy the inequalities 

\\Vt{y,T) - C{t,y,z,p)\\ < \\y - z\\ , 

M 

Utiy.p)-'ip{t,y,z,p)\\<—e \\y-z\\, 

||C(^, 2 /,^\v^) - C(^, 2 /, 2 (^V?)|| < 2 e"^'^* 11 ^^ - ^^11 , 

||^(^, 2 /,^\v^) -' 0 (^, 2 /, 2 ;^</ 2 )|| < 


M 


,-£7i |U1 




and the equality C( 0 , y, z, p) = z. 

To make further notations shorter where it does not lead to confusion, for functions of t, y, z, p 
we will only indicate dependency on the time variable t and write rjt, (ft, Q and fjt instead of 
rit{y,p), (j)t{y,p), Cit^y^^yP) cind 'if{t,y, z,p) respectively. Besides, Remark [2] allows us to 
assume that N = 5. 

Let us introduce a structure of a complete metric space in Tie having dehned the distance 
between a couple of elements {C, t/'*) £ Tie, i = 1,2 with the formula 
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On DJts we declare mappings 

OO OO 

■= (l>t + j [u]{r],,e) - Cd{Cs,e)]<ls + j e:) - if(C, e:)] ds, 

t t 

t 

:=0*z + £ j nl [F{Cs,e) + e^^^G{Cs,^s,e) - {F^{r]s,e) + e^^^Gy{r]s,(j)sG)) Cs] ds. 

0 

For a fixed bundle {y, z, ip) E x T" an element of the space DJl^ generates a solution of 
system fl22l) if and only if 

Ct = G[C,'4^]{t), ^l^t = J^[C,'4^]{t) Vf>0. (27) 

Indeed, the given element of the space DJls generates a solution if and only if for f > 0 the 
inequalities hold: 

Ct = G[C,i^]{t), 

t 

'Ipt - Gt = 'Go - ^ + j [cziCs,£)-i^{Vs,£)+£^^‘^ {H{Cs,'Gs,£) - H{ys,Gs,£))] ds. 

0 

The second one is obvious, whereas to derive the first one, it suffices to write down the solution 
y = Ct of the linear non-homogeneous system y = eP{t; y, (p, e)y + fit), where 

fit) := £ [F{Ct, £) + e^i‘^G{Ctt 'Gt, £) — + ^^^'^G'yir\t, Gt, £)) Ct] , 

which takes the value z at t = 0. Since Gt — Gt ^ f^i t ^ oo, we obtain the only possible initial 
value 

OO 

Go = p>-J [oj{Cs,£) - i^iVsG) + {H{Cs,Gs,£) - H{ys,Gs,£))](^s, (28) 

0 

which leads to Gt = ^[C^G]{t)- Vice versa, if equalities fl27j) are true, then it is evident that 
t ^ {Ct, Gt) is a solution of system ([22]). 

Let us show that the choose suitable Eq for the mapping 

me3{C,G)^{G,J^) 

to be a contraction. In the following we assume that 2^0 < cr. Suppose that {C,G) ^ Then 
inequality (I2TD yields 

OO 

\\HC,G]{t) - Gt\\ < K j [WCs - Vs\\ + WGs - Gs\\] ds < 

t 

^K[2 \\y - z|| + M£V 2 11 ^ _ ^ 11 ] Q-eit ^ ^ 

~ £7 “ £ 

under condition that 

£0^^ < 2/M. (29) 

If we take into account the fact that 

t 

F = ^ly + £ J [F{Vs, b) + £^^‘^G{ys, Gs, b) - {FyiVs, b) + e^^Gyi^s, Gs, b)) ys] ds 


+ 


Me, 


3/2 


e-"^* \\y - ;7|| < —\\y - zl 


M 

£ 


0 
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and inequality fl25l) with = (s, = Vs, = 4^s, h will give us 

ll^[C,^]W - hill < lb - ^11 + 

t 

j [(IICs - Vs\\ + - 0 , 11 ) lie, - Vs\\ + W'lps - 0 .||] ds < 

0 

< e- 2 "^' \\y - ;2|| + \\y - z|| +-—e"^^* \\y - ;2|| < 


7 


< 


1 + 


Men 


+ 


e-"^* \\y - z|| < 2 e-"^* \\y - zj 


as long as 




e 4 

Next, having set Q := ((t, y, z\ (p), 'ipl := 'ip{t,y, z\ip), i = 1,2 and having made the as¬ 
signments in inequalities (1241) 1 and (1251) 1 {y\p>’‘) = * = 1 ) 2 , = rjs, = 4>s, we 

get 

oo 

< f [|lci - CJII + IIU - lAfll] ds < 


< [<='’” IlCi - CflO + [e-- llV-i - .ifll] < 

s>0 ^ s>0 


< 


1 Me 


3/2 


M 




2 4 

3^"* [\\G[C\i^^m - ^[C^b^]b)||] < Ik' - + 


+e-^'>^eK 


^e-ys 


[lie -hs|| + 2||C-r/,|| 114^-0.11] 


sup [e^^^||C]-C 2 ||] + 


s>0 


^e-S7i^5/2^ 


e"^Ms 


+eK {6e + e'^^bsup [e"^" ||Ci 


sup [e^ 7 *|| 0 i- 02 ||] < 11^1 
Ke^/^ 


s>0 


s>0 


+ 


7 s>o 


sup [e^^^|| 0 i- 02 ||] < 


< 




3Meo 

1 +-2 + o +—7^ 

e 2e 4 

If we choose £o > 0 that satishes inequalities (129|) and 

3M£o , Mhl'^ Mhl'^ 

e ^ 2e ^ 4 - ’ ^ ^ 

and also include that the improper integral in the expression for IF converges uniformly in 
y, z, if, then for all e G (0, £o) and {(, 0) G 971, we will have 

WI, 3 (C,i/>) (SIC.V'l.UC.V']) e 


(31) 
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Next, having defined for every ^ ‘^e fhe functions Q = Cit,y, z,(p), ijjl = 

z, if) and having set in inequalities fl 2 ^ and fl25|) (y\(p^) = (C*,-!/)*), i = 1,2, y^ = rjs, 
(p^ = (ps, we draw in entirely the same manner as before that 


^e'yt 


[||UctV'‘l(«)-UctUl(«)||] < 


. 3/2 


.dj. 


< ^sup [e'3- lie] - C]||] + [e'"' ||U - t 

tie ^ tie 

e'’‘[||a[c‘,V'‘iw-ftlctV'"i(()||] < 

< eK (6£ + e-'3‘isup [e”* ||e] - e]||] +-sup [e'3< ||v.] - V]||] < 

d tie 


Ue 


< 


Msq 

4e 


6 + sup [e^^* ||Cl “ Cs ||] + 


Mee 


1/2 


U, 


sup [e'3> llv] - 


u, 


It yields 


/Vfp3/2 , /\y2pl/2 

< y ° »up [e'3- IIU - ■e.]||] + 

tie 


1 /2^ 


EqM ( 6 + Meq 

4+ S 


Msup [e'3‘ lie - C 


u, 


One can easily ensure, that if Eq is small enough, inequalities fl2^ . (150]) grant the validity of 
the contraction conditions in for mapping (13ip for all e G (0,£o)- The fixed point (C,V’) of 
this mapping is precisely the solution of system (| 22 |) which belongs to the space 

By (l28p we can explicitly determine 

0{y,z,p:) := ^{D,y,z,p:). 

Clearly, the mapping d{-,-,-) is continuous in all variables and is Lipschitz-continuous in with 
the constant M/e. Moreover, since C,{ll,y, z,ip) = z, we have 

C(C y, z, (p) = yt {z, e{y, z, (p)), ij{t, y, z, (p) = 0 * {z, e{y, z, (p)). 

Finally, it follows from the very construction of the mapping 'd^(-) := 6{y,z,-) that for all 
{y, z, (p) E Be X T” the equalityo ’dy{(p) = ip holds. Indeed, by (1261) to the point {z, y, ip') G 
Bs X ¥"■, where ip' = ^/^{ip), there is the corresponding point d{z,y,ip') = ip. Thus, there is 
the continuous inverse mapping ^ (•) = "&/,{■) defined on the set dy (T'^). But then dy (T'^) 
being an open-closed subset of a torus, which itself is an open-closed set, coincides with T". 
Hence, 'dy(-) is a homeomorphism of the torus onto itself. □ 

Since detFy{y*,0) ^ 0, the implicit function theorem says that for sufficiently small £o > 
0 , there exists a unique smooth mapping //*(•): [ 0 ,£o] —t such that i/*( 0 ) = y* and 

F{y^{E),E) = 0 for all £ G [0,£o]- On the ground of invariant tori perturbation theory [7119] ITT] 
nans] let us prove the next proposition. 

Proposition 9. There exists such eo > 0 that for all e G (0,eo) system (12^) has an invariant 
torus Te given by equation y = y^{E) + E^s{ip), where the mapping ^^{■): —)■ satisfies 

the Lipschitz condition with the Lipschitz constant L{e) such that L{e) —)■ 0, and g{E) —)■ 0 when 
£ —?• 0. This torus is a local attractor and it attracts all forward trajectories which start from 
an e:(l — q{e)) -neighborhood of the point y^{E). 
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Proof. Let us switch to a new variable ^ in system fl2^ with the substitution y = ?/*(£) + e^. 
It will give us 

^ + ^/eG{y^{e) + v?,e)] , 

(p = uiyPe) + e) + e^^^H{yPe) + e^, cp, e). 

This system can be written as 

? = + (32) 

ip = a(y,{£),£) + £0{t f, £), 

where 

-{p, = [^y{y*{£)G) - Fy{y*, 0)] ^ + y/eG{yPe) + e^, if, e), 

1 

^ / ^'y ^ ^ e^/‘^H{yPe) + ip, e). 

0 

One can easily make sure that Lemma 2.1 from j7] is applicable to system fl32|) . According 
to this lemma, there exists £o > 0 such that for all e G (0, £o) fhe system has an invariant 
torus dehned by equation ^ = ^eip)) where the mapping ^e(-) possesses all of the properties 
mentioned earlier. 

Hence, system 0221 ) in an e^(e:)-neighborhood of point yPe) has an invariant torus 7^ given 
by equation y = y^e) + ef^e{.p)- Let us assign Ze{(p) := y*{£) + ef^e{.p)- For sufficiently small £o 
and e G (0,£o) the point y*{e) lies in together with its ^-neighborhood and for any 

such I/O that ||i/o — 2/*(^)|| < ^(1 ~ ^?(^)) ^md arbitrary ip the equalities 

\\yo - Ze{p)\\ < \\yo - y*{£)\\ + lb*(e) - Ze{p)\\ < £(1 - ^(e)) + eQ{e) = e. 

are valid. Thus, if yo is an arbitrary point of an £(1 — ^/(e))-neighborhood of the point y^e), 
then (i/o, Ze{p)) G Be for all p G ¥"■. Proposition |8] implies that for alH > 0 we have 

ht{yo,Po) -r]t{ze{p),9{yo,Ze{p),Po))\\ < \\yQ - Ze{p)\\, 

1101 {yo,Po) - (t)t{ze{p),e{yo,Ze{p),Po))\\ < \\yo - Ze{p)\\ ■ 

It means that the forward trajectory of the point (i/o, po) is attracted to the forward trajectory 
of the point {ze{p),6{iiQ, Ze{p),po)). To ensure that the latter trajectory lies on the invariant 
torus, the hxed point condition has to be fulhlled; 

e{yo,Ze{p),Po) = p. 

Let us show that such a hxed point on the torus exists. The number Eq may be taken small 
enough for ML{e) < 1 to hold for all £ G (0,eo)- Then, in accordance with Proposition [HI for 
arbitrary pi,p 2 G T” we have 

\\e{yo,Ze{pi),Po) - e{yo,Ze{p2),Po)\\ < ML{e) \\pi - p2 \\. 

Consequently, the contraction mappings principle implies the existence of the unique point 
= p*{yo, Po) G T” such that 6*(i/o, Ze{p*), Po) = P*, which means that the forward trajectory 
of the point (yo,Po) is attracted to the forward trajectory of the point (ze(pP, p^,) of the 
invariant torus %. At the same time 

\\yt{yo,Po) - Ze (0i(</5*))|| < 2e"^^‘ \\yo - Ze{pP)\\ , t > 0. (33) 

Let us note that since can be found using the method of subsequent approxima¬ 
tions, pPyQ^po) continuously depends on {iio,Po). It follows from (12^ for f = 0 that 
9{ze{po),Ze{pQ),Po) = Po, and therefore, p^{ze{po),Po) = Po- Thus, pPyo,po) Po when 
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I/O —>■ Ze{(po)- But then at the same time Zs{(p^:{yo, cpo)) —)■ yo. By fl3^ this yields that for 
any A > 0 there exists such 5 > 0, that the forward trajectory of a point which lies in a 
5-neighborhood of the torus 7^ belongs to a A-neighborhood of this torus and is attracted to 
the latter. This does mean that 7^ is a local attractor. □ 

Remark 3. In case of a quasi-periodic flow on an invariant torus, an estimate similar to fl33|) 
was obtained in na. 

Proposition 10. The attraction basin of the invariant torus 7^ contains the forward invariant 
set ([0, \\ne^\]) n 5J x x T*". 

Proof. Turning back to Propositions [6] and [71 it is enough to prove that Bf'^^{y*) x T” lies in 
the attraction basin of the torus 7^. Suppose that (i/o, (po) is an arbitrary point of the domain 
B'f'^'^{y*) X T”. Let us choose a hnite series of points {yi]\^i such that ||i/i-i — i/i|| < e with 
the last point yj lying in an e:(l — p(e))-neighborhood of the point y^{e). Then, applying 
Proposition [8] step by step, we can prove that there exists such (fj e that the forward 
trajectory of the point ((fio,yo) is attracted to the forward trajectory of the point {ipi,yi), 
which in turn, according to Proposition [9| is attracted to the torus %. Thus, under action of 
the semi-flow of system (|22|) {(po,yo) is attracted to the torus %. □ 

7. Conclusion 

In this paper, we have analyzed a kind of transient processes that are observable in a fast- 
slow system in a neighborhood of an i. m. s. m. and that can be interpreted as a dynamical 
bifurcation of multi-frequency oscillations. The change in the phase variables x{f) behavior 
— the switch from damping oscillations to the multi-frequency ones, which are asymptotically 
close to motions on the invariant torus % — is caused by the slow evolution of the parameters 
u{t), which results in transition of the latter from the stability zone Bg to the zone of complete 
instability Bu- 

We should note that there is certain connection between the obtained results and the the¬ 
ory of bifurcations without parameters (see [29] and references there). In [29] system ([3]) was 
investigated in the case where the i. m. s. m. a; = 0 consists completely of equilibria and the 
spectrum of the operator f^[0,0,e) lies on the imaginary axis. In this situation the transfor¬ 
mations of the phase portrait are caused by different kinds of hyperbolicity of the system of 
the hrst approximation with respect to points of the i. m. s. m. ( 0 ,m) if m 7 ^ 0 . 
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